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Employing Molecular Dynamics simulations of a chemically realistic model of 1,4-polybutadiene 
between graphite walls we show that the mass exchange between layers close to the walls is a slow 
process already in the melt state. For the glass transition of confined polymers this process competes 
with the slowing down due to packing effects and intramolecular rotation barriers. 
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Polymer dynamics in confinement is both of funda- 
mental interest concerning our understanding of the glass 
transition 0, 0| , as well as of high technological impor- 
tance for the performance of composite materials Q. In 
experiments 0-Q| as well as simulations l8l-flfj] varying 
effects of confinement on the glass transition and even 
conflicting results have been observed for one and the 
same polymer- wall system |J-|6|. Obviously, the dynam- 
ics will be inhomogeneous as well as anisotropic in the 
vicinity of the confining surface, but it is not clear how 
large this effect is and how its extent into the bulk of the 
polymer is connected with structural anomalies close to 
the surface. These in turn will sensitively depend on the 
details of the polymer-solid interactions. Thus we found 
it highly desirable to perform a detailed simulation anal- 
ysis of a well defined chemically realistic polymer-wall 
system and take advantage of the strength of the simu- 
lation approach to provide insight not (yet) obtainable 
from experiment. The segment length of the polymer 
and the radius of gyration of the chains have been iden- 
tified as the relevant length scales over which the wall 
effect changes different aspects of the polymer dynamics 
0, i.e., length scales from the sub-nanometer scale to at 
most 10 nm. When one tries to model the mechanical 
properties of polymer composites on the other hand, one 
often has to assume a region of modified viscoelasticity 
within the polymer of several hundred nanometers ex- 
tent H . The origin of the discrepancy between these two 
findings is not understood. 

In bulk glass forming polymers one has a competition 
between two mechanisms leading to a time scale sepa- 
ration between vibrational and relaxational degrees of 
freedom, the caging effect as captured by mode coupling 
theory which is essentially a hard sphere packing 
phenomenon, and the time scale separation introduce d by 
the presence of rotational barriers in the chains [lU, Il3j . 
The latter mechanism is absent in the widely used coarse- 
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grained models (J Q, Il5j . Both mechanisms compete 
upon approaching the glass transition in polymer melts. 
As we will show, attractive confining walls introduce an- 
other mechanism for time scale separation. This mecha- 
nism is the slow desorption kinetics of monomers at the 
surface (which in our case is graphite) leading to a slow 
layer exchange dynamics on the scale of (at least) the 
radius of gyration of the chains. 

Simulation method We have performed Molecular 
Dynamics (MD) simulations of a well- validated, chem- 
ically realistic model |16l4l9f of a 1,4-polybutadiene 
(PBD) melt at a graphite surface. Unlike other studies 
using polyethylene [20l l2~lj |. a polymer that easily crys- 
tallizes, we focus here on a glassforming polymer which 
has been extensively characterized in the bulk. The cur- 
rent model quantitatively reproduces all available exper- 
imental results on PB bulk dynamics (l6l - [l9| and its 
glass transition temperature of about 180 K. The (0, 0, 1) 
graphite surfaces are represented by layers positioned at 
z = 0, z = D = 9.97 nm, and z = 10.31 nm with a pe- 
riodic boundary condition imposed at L z = 10.65 nm. 
All graphite carbon atoms are kept at fixed positions 
in the simulation. The force field for graphite is taken 
from and Lorentz-Berthelot combining rules for the 
Lcnnard-Jones interactions are used. In the x and y di- 
rections periodic boundary conditions are employed at 
L x = 14.91 nm and L y = 14.76 nm, respectively. The 
density at the center of the simulation box is adjusted 
close to its bulk value at the respective temperature, 
and the melt structure is equilibrated with a sequence of 
Langevin dynamics simulations at reduced excluded vol- 
ume under NVT conditions and MD post equilibration. 
Finally, we performed NVT MD using the GROMACS 
package (23J. 

Results The dispersive interactions at the interface 
between both materials lead to a strong layering of the 
polymer segments at the graphite surface. Furthermore, 
for the chains of 29 repeat units which we have simulated, 
there is also an observable layering on the scale of the 
whole polymer chain due to the correlation hole effect 
24] . In Fig. T the monomer density and the density of 
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chain center of mass positions are shown for T = 353K 
(about 2T ), normalized to their bulk values. 




z, nra 



FIG. 1: Layering effects in the density (dashed line) and 
chain center of mass density (connected circles) of a confined 
polybutadiene melt. Both curves are normalized by the cor- 
responding densities in the bulk. 

We checked that the density dependence we obtained 
in the simulation was symmetric with respect to the dis- 
tance to the two walls to judge the quality of equilibra- 
tion. Fig. Q] displays half of the symmetrized profile. The 
length scale involved in the density layering is the size of 
the monomers (width of the chains), a ~ 0.5 nm. The 
perturbation of the melt structure propagates for about 
2.5mn from the graphite surface. Over the same distance, 
we observe a layering effect in the center of mass positions 
of the chains. The maximum in this density lies around 
1 nm which is slightly smaller than the radius of gyration 
of the chains in the bulk because chains within this region 
of increased density are oriented parallel to the surface, 
i.e., their gyration tensor ellipsoid orients with its long 
axes in the x, y-plane. In contrast to the bulk behavior, 
at a wall the correlation hole effect therefore introduces 
a well defined liquid like layering in the center of mass 
density. 

Both of these layering effects are connected with in- 
homogeneous dynamics in the film which we will trace 
to the slow adsorption/desorption process of monomers. 
Let us first focus on the orientational dynamics on the 
segmental scale as would be observable in nuclear mag- 
netic resonance (NMR) experiments. We will study the 
following correlation function [25[ 

C(t) = (P 2 (cos[0(i)]) P 2 (cos[0(O)])) , (1) 

where 8(t) is the angle between the orientation of the dou- 
ble bond in the polybutadiene monomer and a given mag- 
netic field orientation, P 2 is the second Legendre polyno- 
mial and the angular brackets indicate a thermal average. 
We have determined this correlation function in a layer 
resolved manner with a total of 8 layers between the two 
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FIG. 2: Inhomogeneous rotational dynamics of bonds as 
showing up in an NMR correlation function C(t) defined in 
the text. We show the relaxation function for two different 
layers in the film and the average relaxation in the film for 
an orientation of the external field parallel to the walls (X) . 
We compare this to the bulk relaxation and the average film 
relaxation for an external field oriented perpendicular to the 
walls (Z). 



graphite surfaces at z — and z = D. A given double 
bond is contributing to the relaxation function for layer 
j if its center is within this layer at t = 0. Most of the 
relaxation functions are shown for an orientation of an 
external magnetic field parallel to the graphite surface 
(x-direction) . One can clearly observe that the dynam- 
ics is slowed down strongly, going from the center of the 
film (blue, dash-dotted line) to layer 1 (green, dashed 
line) directly at the walls. The relaxation in the center 
of the film is almost the same as in bulk (blue, dotted 
line) which was determined separately. The average sig- 
nal from the thin film is shown by the black dashed curve 
for a field orientation parallel to the walls and by the 
black continuous curve for a field orientation perpendic- 
ular to the walls (z-direction). For the x-orientation the 
relaxation of the last percent of correlations is about an 
order of magnitude slower than in the bulk, whereas for 
the z-direction this relaxation appears to be slowed down 
by 2 to 3 orders of magnitude! For both relaxation func- 
tions a final slow process seems to be setting in around 
2 ns. 

In the bulk such a process (around 2% of the decay) 
in the C — H bond reorientation was shown to be due 
to a coupling of local to large scale relaxation, i.e., con- 
formational relaxation of the whole chain [20]. At the 
wall we are here observing a similar coupling which is 
furthermore most prominent for the perpendicular reori- 
entation. To quantify the large scale motion of the chain 
in the vicinity of the wall we are presenting in Fig. 3 the 
mean square displacement of the center of mass of the 
chains distinguished with respect to directions parallel 
and normal to the walls. This time we have resolved 7 
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FIG. 3: Mean-square-displacement of the center of mass of 
chains for different layers chosen to have a thickness equal to 
the bulk R g . We differentiate between the average displace- 
ment parallel to the walls, xy, and perpendicular to the walls, 
z. 



layers, i.e., every layer has a thickness of the order of 
the bulk R g ~ 1.4 nm. A chain contributes to the dis- 
placement function of layer i when its center of mass is 
in that layer at t — 0. In the bulk, all displacements in 
different Cartesian directions are equal and the same is 
true for the displacement we observe in the middle of the 
film (layer 4, blue and red continuous lines). However, 
due to a slightly larger density than for the bulk simu- 
lations, the displacement in the center of the film is still 
about 30 % slower than in the bulk. We can quantify 
this by using an estimate of the Rouse time, defined by 
{{AX cmtl (T R )f) = R 2 g /3 ~ 0.7 nm 2 . For the bulk this 
yields rjj ulk ~ 9 ns and for the center of the film we obtain 
T middic ^ yi ns. At the wall, in layer 1, the motion of the 
chain is not yet isotropic on a displacement scale of R g . 
For the displacement parallel to the graphite surfaces we 



read off a Rouse time ~ 15 ns, whereas for the dis- 
placement perpendicular to the walls it is rjj' z ~ 39 ns. It 
is interesting to note that from about 1 ps to about 50 ps 
the displacement parallel to the graphite walls in layer 
1 is actually slightly faster than the corresponding dis- 
placement in the center of the film. Torsional transitions, 
which are the relevant local motions on these time scales 
are, actually slowed down close to the wall (not shown), 
so we take this as an indication that dynamical correla- 
tions between local motions are changed in the vicinity 
of the wall. 

One could be tempted to describe the difference be- 
tween parallel and perpendicular mobility by assigning 
different friction constants to the motions in different di- 
rections, however, such an approach fails to account for 
the short time (1 to 10 ps) plateau of the perpendicular 
displacement. This plateau signifies a time scale separa- 
tion between short time vibrational motion (up to 1 ps) 
and longer time relaxational motion. 

In the bulk, caging due to packing effects and in- 



tramolecular barriers sets in around T = 273 K [12J 
and both effects can not be the reason for the plateau 
in the perpendicular displacement, as we have observed 
that the parallel motion close to the wall is as fast as 
the motion in the bulk. The new slow process occurring 
with polymers at attractive walls is the adsorption and 
desorption kinetics [13, HH of monomers at the wall, and 
phenomenologically it gives rise to the same two-step re- 
laxation behavior as the other two mechanisms. 
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FIG. 4: Adsorption autocorrelation function for single 
monomers. The assignment to layers is done on the basis 
of the position of the center of mass of the chain at t = 0, 
where by definition at t = the considered monomer has to 
be attached to the graphite surface. 



We quantify the desorption kinetics by the following 
correlation function 



(s(t) S (0)) - ( S ) 2 _ (s(t)) - (a)* 

( s (o>(o)> - ( S ) 2 i - <*> 2 



(2) 



where s(t) = 1 if the monomer is adsorbed (i.e., is within 
the first layer observable in the density profile) at time 
t and s(t) = otherwise, and where we have employed 
the constraint s(0) = 1 for the collection of the autocor- 
relation function. In Fig. |4] we can see that this auto- 
correlation function decays on a time scale comparable 
to the relaxation time defined from the mean squared 
center of mass displacement in z-direction (see Fig. [3]). 
Furthermore, especially for the monomers at the wall, a 
clear two-step decay is observable with a plateau extend- 
ing over the same time interval as the one observed for 
the center of mass msd and a plateau height above 0.95. 
We also observe that this plateau extends further at the 
lower temperature of T = 323 K which we included in this 
figure. Defining a time scale for the desorption kinetics 
by <P(t w ) — 0.5, we read off t w = 10 ns at 353 K and 
t w = 38 ns at 323 K. Their ratio of 3.8 is larger than the 
ratio of 2.4 we found between the structural relaxation 
times in the center of the film for these two tempera- 
tures, which in turn agree with the corresponding bulk 
relaxation times [l9[ at the same density. Therefore, al- 
ready from our simulations well above T g we obtain an 
indication that the presence of this additonal slow pro- 
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cess may have a strong influence on the glass transition 
of this confined system. 

Conclusions Confining 1,4-polybutadiene between 
graphite walls creates layering effects in the density as 
well as the center of mass density extending to about 2Rq 
away from the walls. Both structural features influence 
the way in which different correlation functions depend 
on the distance to the walls as well as the way in which 
they become anisotropic. We have shown for a double- 
quantum NMR signal that this anisotropy as well as the 
slowing down induced by the walls become observable ex- 
perimentally in the final few per-cent of the decay. The 
anisotropy can be very large and should be observable in 
experiments of strongly confined systems where the di- 
rection of the probing field with respect to the walls can 
be controlled. In other relaxation functions, e.g., scatter- 
ing experiments, averaging over different directions of the 
momentum transfer and over the whole confined sample, 
might lead to only small residual observable effects, al- 
though the alteration of the dynamics in the vicinity of 
the wall can be rather large. Experimentally it is there- 
fore highly desirable to develop techniques which can re- 
solve relaxation behavior as a function of distance to the 
walls with a resolution on the nano-meter scale and to 
perform orientationally resolved experiments. Our sim- 



ulations should be a guide to judge what is needed in 
this context, and hence allow to clarify the existing con- 
troversies of interpretation. We have shown that the 
slow desorption kinetics constitutes a third mechanism 
for time scale separation in glass forming polymers in 
addition to packing effects and conformational barriers 
present in bulk samples. How the competition between 
these mechanisms develops upon approaching the glass 
transition in different samples is a very interesting ques- 
tion which deserves further experiment and simulation 
studies. One might speculate that the extend over which 
the walls influence the visoelastic response of a confined 
melt Q increases with decreasing temperature and in- 
creasing chain length (R g dependent layering). However, 
to address such a question in simulations, a mapping of 
our atomistic model to a coarse-grained model which can 
then be studied on larger length scales needs to be per- 
formed, which remains as a challenge for future work. 

Acknowledgment: We acknowledge funding through 
the German Science Foundation through the focused 
funding program SPP 1369, sub-project C2. We are 
grateful to the Jiilich Supercomputer Center for com- 
puter time on the JUGENE and JUROPA computers 
through project HMZ03, and to the European network 
of excellence SoftComp for computer time. 



J. Baschnagel and F. Varnik, J. Phys.: Condens. Matter [15 
17, R851 (2005). 

M. Alcoutlabi and G. B. McKenna, J. Phys.: Condens. 
Matter 17, R461 (2005). [16 
S. G. Advani, Processing and Properties Nanocompos- 
ites (Worlds Scientific, Singapore, 2006); D. Gay and S. [17 
V. Hoa, Composite Materials, Design and Applications 
(CRC Press, Boca Raton, 2007) [18 
A. N. Raegen, M. V. Massa, J. A. Forrest and K. Dalnoki- 
Veress, Eur. Phys. J. E 27, 375 (2008). [19 
A. Serghei, H. Huth, C. Schick and F. Kremer, Macro- 
molecules 41, 3636 (2008). [20 
S. Kim, S.A. Hewlett, C.B. Roth and J.M. Torkelson, 
Eur. Phys. J. E 30, 83 (2009). 

P. A. O'Connell and G. B. McKenna, Science 307, 1760 
(2005). [21 
G .D. Smith, D. Bedrov and O. Borodin, Phys. Rev. Lett. 
90, 226103 (2003) [22 
J. A. Torres, P. F. Nealey and J. J. de Pablo, Phys. Rev. [23 
Lett. 85, 3221 (2000) 

H. Eslami and F. Miiller-Plathe, J. Phys. Chem. B 114, [24 
387 (2010). 

W. Goetze, Complex Dynamics of Glass-Forming Liq- [25 
uids: A Mode- Coupling Theory, (Oxford University [26 
Press, Oxford, 2008). 

S. Krushev, W. Paul, Phys Rev E 67, 021806 (2003). [27 
M. Bernabei, A. J. Moreno and J. Colmenero, Phys. Rev. 
Lett. 101, 255701 (2008). [28 
[14] T. Pakula, J. Chem. Phys. 95, 4675 (1991). 



J. Baschnagel, K. Binder, M. Milchev, in Polymer Sur- 
faces, Interfaces and Thin Films, A. Karim and S. Kumar 
Eds., (World Scientific, Singapore, 2000). 
G. D. Smith and W. Paul, J. Phys. Chem. A 102, 1200- 
1208 (1998). 

G. D. Smith, W. Paul, M. Monkenbusch, and D. Richter, 
J. Chem. Phys. 114, 4285-4288 (2001). 
W. Paul, D. Bedrov and G. D. Smith, Phys. Rev. E 74, 
021501, (2006). 

W. Paul and G. D. Smith, Rep. Prog. Phys. 67, 1117- 
1185 (2004). 

K. C. Daoulas, V. A. Harmandaris and V. G. Mavrantzas, 
Macromolecules 38, 5780 (2005). V. A. Harmandaris, K. 
C. Daoulas and V. G. Mavrantzas, Macromolecules 38, 
5796 (2005). 

O. Alexiadis, V. Mavrantzas, R. Khare, J. Beckers and 

A. Baljon, Macromolecules, 41, 987 (2008). 
W. A. Steele, Surf. Sci. 36, 317 (1973). 

B. Hess et al. J. Chem. T heory Comput. 4, 435 (2008); 
|http://www . gromacs.org 

P. G. De Gennes, Scaling Concepts in Polymer Physics, 
(Cornell University Press, Ithaca, 1979). 
K. Saalwachter, Prog. Nucl. Mag. Res. Sp. 51, 1 (2007). 
W. Paul, G. D. Smith and D. Y. Yoon, Macromolecules 
30, 7772 (1997). 

K. A. Smith, M. Vladkov, and J.-L. Barrat, Macro- 
molecules, 38, 571 (2005). 
M. Erber et al., Polymer 51, 129 (2010). 



